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GLOBAL EXISTENCE AND SEMICLASSICAL LIMIT FOR QUANTUM 
HYDRODYNAMIC EQUATIONS WITH VISCOSITY AND HEAT 

CONDUCTION 

XUEKE PU AND BOLING GUO 


Abstract. The hydrodynamic equations with quantum effects are studied in this paper. 
Eirst we establish the global existence of smooth solutions with small initial data and 
then in the second part, we establish the convergence of the solutions of the quantum 
hydrodynamic equations to those of the classical hydrodynamic equations. The energy 
equation is considered in this paper, which added new difficulties to the energy estimates, 
especially to the selection of the appropriate Sobolev spaces. 


1. Introduction 

The hydrodynamic equations and related models with quantum effects are extensively 
studied in recent two decades. In these models, the quantum effects is included into the 
classical hydrodynamic equations by incorporating the first quantum corrections of 0{h'^), 
where h is the Planck constant. One of the main applications of the quantum hydrody¬ 
namic equations is as a simplified but not a simplistic approach for quantum plasmas. In 
particular, the nonlinear aspects of quantum plasmas of quantum plasmas are much more 
accessible using a fluid description, in comparasion with kinetic theory. One may see the 
recent monograph of Haas [8] for many physics backgrounds and mathematical derivation 
of many interesting models. Many other applications of the quantum hydrodynamic equa¬ 
tions consisting of analyzing the flow the electrons in quantum semiconductor devices in 
nano-size [7], where quantum effects like particle tunnelling through potential barriers and 
built-up in quantum wells, can not be simulated by classical hydrodynamic model. Similar 
macroscopic quantum models are also used in many other physical fields such as superfluid 
and superconductivity [6]. 

Let us first consider the following classical hydrodynamic equations in conservation form, 
describing the motion of the electrons in plasmas by omitting the electric potential 


( dn dn, 

dt m dxi ’ 

(I.la) 


(1-lb) 

dW d 

^ Qx ^j^ij ^i^ “0; 

(1.1c) 


where n is the density, m is the effective electron mass, u = (ui,U 2 ,M 3 ) is the velocity, Hj 
is the momentum density, Pij is the stress tensor, W is the energy density and q is the 
heat flux. In this system, repeated indices are summed over under the Einstein convention. 
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This system also emerges from descriptions of the motion of the electrons in semiconductor 
devices, with the electrical potential and the relaxation omitted. 

As in the classical hydrodynamic equations, the quantum conservation laws have the same 
form as their classical counterparts. However, to close the moment expansion at the third 
order, we define the above quantities H^, and W in terms of the density n, the velocity u 
and the temperature T. As usual, the heat flux is assumed to obey the Fourier law q = —kVT 
and the momentum density is defined by H^ = mnui^ where m is the electron mass and u 
the velocity. The symmetric stress tensor P^- and the energy density W are defined, with 
quantum corrections, by 


Pij — JlTdij 


h?n 

12 m dxidx. 


■ logn + 


and 

W = ^nT + \mn\u\^ — logn + 

2 2 24m 

respectively, where h is the Planck constant, and is very small compared to macro quantities. 

As far as the quantum corrections are concerned, the quantum correction to the energy 
density was first derived by Wigner |24) for thermodynamic equilibrium, and the quantum 
correction to the stress tensor was proposed by Ancona and Tiersten [5] and Ancona and 
lafrate [1] on the Wigner formalism. See also [7] for derivation of the system dm by a 
moment expansion of the Wigner-Boltzmann equation and an expansion of the thermal 
equilibrium Wigner distribution function to leading to the expression for H and W 

above. We also remark the quantum correction term is closely related to the quantum Bohm 
potential [4] 


Q(n) = - 


Ay/n 
2 m ^Jn ’ 


where n is the charge density. It relates to the quantum correction term in P^- with 

^2 ^2 ^2 

—nVQ{n) = - —div(nV^logn) = -—AVp -div(V-v/ri ® a/u)- 

4m 4m m 

For the system inj, there is no dissipation in the second equation. Given n and T, the 
second equation if hyperbolic, and generally we can not expect global smooth solutions for 
this system. In this paper, we consider the following viscous system by taking into account 
the stress tensor §, 


' dtn + V ■ (nu) =0, ( 1 . 2 a) 

1 1 

dtu + u ■ Vu H- V(nT) — — —^div{n(V G V) logn} =-divS, (l-2b) 

^ mn Izm^n mn 

dtT + u-VT+ §PV • u - J-V • (kVT) + 7 -^—V • inAu) = :r^{V • (uS) - u ■ divSj. (1.2c) 

V O oTl oOTTZTZ oTRTl 

Here, S is the stress tensor defined by 


§ = piVu + iVu)^) + A(divM)I, 


where I is the dx d identity matrix, /i > 0 and A are the primary coefficients of viscosity and 
the second coefficients of viscosity, respectively, satisfying 2/r + 3A > 0. Without quantum 
corrections (i.e., setting /i = 0 ), this system is exactly the classical hydrodynamic equations 
studied in the seminal paper of Matsumura and Nishida m- 

Although important, there is little result on the system (EH) to the best of our knowl¬ 
edge. But there does exist a large amount of work for system very similar to (11.21) . These 
work comes from two main origins. The first one is from the quantum correction to various 
hydrodynamic equations, especially in semiconductors and in plasmas. Gardner [7] derived 
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the full 3D quantum hydrodynamic model by a moment expansion of the Wigner-Boltzmann 
equation. Hsiao and Li m reviewed the recent progress on well-posedness, stability analysis, 
and small scaling limits for the (bi-polar) quantum hydrodynamic models, where the inter¬ 
ested readers may find many useful references therein. Jungel m proved global existence 
of weak solutions for the isentropic case. See also [TMTB] . 

The other one, being equally important, emerges from the study of the compressible fluid 
models of Korteweg type, which are usually used to describe the motion of compressible 
fluids with capillarity effect of materials. See Korteweg [15] and the pioneering work of 
Dunn and Serrin |5]. The reference list can be very long, and we only mention a few of them. 
Hattori and Li mn] considered the local and global existence of smooth solutions for for the 
fluid model of Korteweg type for small initial data. Wang and Tan [53] studied the optimal 
decay for the compressible fluid model of Korteweg type. Recently, Bian, Yao and Zhu |3] 
studied the global existence of small smooth solutions and the vanishing capillarity limit of 
this model. Jungel et al [12j showed a combined incompressible and vanishing capillarity 
limit in the barotropic compressible Navier-Stokes equations for smooth solutions. 

Almost all of the above mentioned results considered the isothermal case, studying only 
the continuity equation and the momentum equation, or with electric potential described 
by a Poisson equation. To the best of our knowledge, there is no mathematical studies for 
the full quantum hydrodynamic system (HID. The system (USD is itself interesting, since 
the energy equation also includes the quantum effects through the energy density W, which 
brings new features into this system. This makes it different from the previous known results, 
to be precisely stated in the following. 

The aim of this paper is two fold. On one hand, we show the global existence of smooth 
solutions for (HID with fixed constant h > 0 when the initial data is small near the con¬ 
stant stationary solution {n,u,T) = (1,0,1). To be precise, we denote the perturbation by 
{p,u,9) = {n — l,u,T — 1) and transform the (11.21) into (12.11) . The result is then stated in 
Theorem 12.21 for (12.11) , where the estimates is stated in terms of the planck constant h, and 
we can see clearly how the quantum corrections affect the estimates. On the other hand, 
since (USD modifies the classical hydrodynamic equations to a macro-micro level in the sense 
that it incorporates the (micro) quantum corrections, it is expected that as the Planck con¬ 
stant —>■ 0, the solution of the system (USD converges to that of the classical hydrodynamic 
equations. This limit is rigorously studied in this paper and stated in Theorem 12.31 In 
particular, algebraic convergence rate is given in terms of h. 

Among others, one of the main novelties is the selection of the Sobolev space like (p, u, 9) G 
jjk +2 ^ X . The underlying reason lies in the fact that the quantum effects in the 

energy density involves higher order derivatives of the velocity, and hence we cannot seek 
solutions in the same Sobolev spaces for 9 and u. 

This paper is organized as follows. In Sect. 2, we present some preliminaries. We translate 
the system (11.21) into a convenient form and state the main results in this paper. In Sect. 3, 
we give the a priori estimates, and then prove Theorem (12.21) (existence result) at the end of 
this section. Finally, in Sect. 4, we prove Theorem (12.31) . by showing the convergence of the 
solutions of the quantum hydrodynamic equations dUD to that of the classical hydrodynamic 
equations (USD- The algebraic convergence rate is also given in terms of the Planck constant 
h. 

Notations. Throughout, C denotes some generic constant independent of time f > 0 and 
the Planck constant h> Q. Let p G [l,oo], denotes the usual Lebesgue space with norm 
II • \\lp- When p = 2, it is usually write || • || = || • ||lp, omitting the subscript. Let denote 
the Sobolev space of the measurable functions whose generalized derivatives up to order 
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belong to L? with norm || • ||^fc = (Xi=o ‘ denotes the homogeneous Sobolev 

spaces and \A^ B] = AB — BA denotes the commutator of A and B. 

2. Preliminaries and Main Results 

In this section, we reformulate the system in convenient variables. First we take 

(n, u, T) = (1,0,1) to be a constant solution to (11.21) and consider 

p = n—l,u = u,0 = T—1. 

In these unknowns, with m = 1, the & transforms into 
dtp + u ■ Vp + (1 + p)divu = 0, 

dtu - - —Am — ^ VdivM = —u ■ Vm — V0 — ^ ^ Vp 

p+1 P+1 P+1 

AVp div(Vi/p + 1 ® VVp + 1) 

< "^12^ + 1 ~ y p +1 ’ 

2k 2 

with initial data 

(p, M, e){0, x) = {po, Mo, 6 »o)(a;) = (no - 1, mq, Tq - l)(a::). 

We first state the local-in-time existence of smooth solutions to (E+J. To be precise, we first 
set 

|||(p,m,P)||| 2 :=||(p,M,P)||i. + llVpIli. -f ||(Wp,ftVM)||i. + Ilft^Aplli., 

|||(p,M,P)|||^=|||(p,u, 0 )lll^l + lll(vV,vVv^P)|||g, 

and 

£:fc(0,T) = i (p,M,P)(t, •) : sup |||(p,M,P)(t)|||fc < oo 
( te[o.T] 

Theorem 2.1 (Local existence). For any initial data such that uq > 5 > 0 is satisfied and 
(Pq = no — l,uo,6o) S x x (k > 3), there exists some T > 0 sueh that the 

Cauehy problem m has a unique solution {p,u,9) in [0,T] such that {p,u,9) € £k(0,T) 
and 

lll(P:^^,^')(i)lllfc < Ck\\\{po,uo,9o)\\\l. 

This theorem can be proved in a similar fashion as in by the dual argument and 
iteration techniques, and hence omitted for brevity. 

Now, we consider the global existence of smooth solutions. Let T > 0, we set 

E= sup |||(p,M, 6 »)(t,-)||| 3 - (2.2) 

0<t<T 

One of the main purpose is to show the following 
Theorem 2.2 (Global existence). Suppose the initial data 

(p,m,6»)(0) G X H'^ X 

and set Eq := |||(p, m, P)(0)| jja < oo. There exists some Hq > 0, Eq > 0, mq > 0 and 
Cq < oo, sueh that if Eq < Sq and h < ho, then there exists a unique global in time 


(2.1a) 

(2.1b) 
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solution (p, u, 9) £ ^3(0, T) of the Cauchy problem (12.11) for any t £ (0, c»), and the following 
estimates hold 

|||(p,w,0)(t)|||2 + uo / ^||V'=(p,u,0,ftVii,Wp)(s)fds<Co|||(p,w,0)(O)|||2, 

•^0 k=i 

where the constants uq > 0 and Cq > 0 are independent of time t and h. 

Formally, as h^o, ^ tends to the following classical hydrodynamic equations for 
(studied in [S]) 

dtp^ + ■ Vp° + (1 + p°)divu° = 0, (2.3a) 

dt9°-^^^^A0° = -u°-Ve°-^i0° + l)V-u° 

+ ^(1^7) {f + A(divu°)2} . (2.3c) 

The convergence result is stated in the following 

Theorem 2.3 (Semiclassical limit). Let {p^,u^,0^) be the solution of (12.11) and (p°,M°,0°) 
be the solution of (12.31) with the same initial data {pQ,utj,0Q) £ x H'^ x H^. Then for all 
fixed time T £ (0, 00), we have the algebraic convergence 

sup IKp'^ - p°,M^ - - 6»°)||^i < [cse'^i'^/ci] ;i'‘ 

te[o.T] 

and 

sup Up^-p°,u^-u°, 0^-0°)\\%,< [c 2 e‘^^^/ci]h'^, 
te[o.T] 

for some constant positive constants ci and C 2 , independent of h and t. 

The following three lemmas will be frequently used, and hence cited here for reader’s 
convenience. 

Lemma 2.4 (Gagliardo-Nirenberg [lO])- Let p,q,r £ [l,oo] and 0 < i,j < I be integers, 
there exist some generic constants 0 £ [0,1] and C > 0, such that 

where 

N N N 

J--=0{l--) + {l-0){^--). 

When 0=1, I — j ^ N/q. 

Lemma 2.5. Let g{p) and g{p,0) be smooth functions of p and {p,0), respectively, with 
bounded derivatives of any order, and ||p||l= < 1- Then for any integer m > 1, we have 

||V”^5 (p)|Up <C||V>|Ur, ||V™5 (p, 0 )|Up <q|V>,V™p|Up, Vpe [l,oo], 
where C may depend on g and m. In particular, 

(^) lU- <C'II(p,0)II^m,., Vpe [1,00]. 

Proof. This can be proved in a similar fashion as in [23] and [5] making use of the Gagliardo- 
Nirenberg inequality, and hence omitted here for brevity. □ 
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Lemma 2.6 (Kato-Ponce [14)1. Let a be any multi-index with |a| = fc andp € (l,oo). Then 
there exists some constant C > 0 such that 


||5“(/5)||l. < C{\\f\\LP^\\9\\Hs.. + ll/llff.,.3ll5lU-}, 

\mJML. < C{||V/|U.,||g||^._,„, + ||/||^...3ll5lU^4}, 

where f^g €S, the Schwartz class andp 2 ,P 3 S (l,+cx 3 ) such that 

1 _ 1 1 _ 1 1 

P Pi P2 PZ Pi ' 

3. A priori estimates 

In this section, we establish useful a priori estimates of the solutions to m- First of 
all, we let the Planck constant < 1. To simplify the proof slightly, we assume that there 
exists a positive number e <C 1 such that 

E= sup |||(p,u, 6 »)(t )|||3 < e, (3.1) 

tG[0,T] 

which together with Sobolev embedding, implies that 

sup \\{p,u,9),V{p,u,9),V^p,hV^{p,u),h‘^V^p\\L<^ < CE < Ce, (3.2) 

t6[0,T] 

and from m the following 

< C||V • (u(l + p))||lp < CF; < Ce, VI < p < oo, (3.3) 

and 

A-k 

<-:^||A 6 >||ip + ||u • V0 ||lp + ||divM||ip + ?i^||V((l + p)Au)||lp 

3 (3.4) 

+ C'||Vu||^ 2 p < CF; < Ce, Vl<p< 6 . 

In particular, we choose e small enough such that 

sup |1 (p,6I)(<)||loo < 1/2. (3.5) 

te[o.T] 


3.1. Basic estimates. Now, we consider the zeroth order estimates for the system (EH). 
As in m , we set 


s — (1 + 9)/ (1 + p)^/^ — 1, 

and define a function E^{p, u, s) for p,u = (u^, u^) and s by 


- 1 - f) + 2-/^ 


,2 p, 3i?(l 
Mf + Rsp H- 


p)s 


(3.6) 

(3.7) 


The following lemma is proved in |19j . 


Lemma 3.1. There exists constants 0 < p 2 < 1/2 and 0 < Ci < C 2 < oo such that E^ is 
positive definite, i.e., 

p 2 + | y |2 ^ ^2 < ^^^0 < ^^(^2 ^ |^|2 ^ 02 )^ 1^1 < 


We first prove the zeroth order estimates in the following 


Proposition 3.1. There exists a constant Eq > 0 such that if E < e < Eq, then the following 
a priori estimates holds for all t G [0,r], 

|||(p,u,6»)(t)|||^ + j/Q [ \\D{p,u,9){s),{hAu,nAp){s)\\'^ds < C\\\{p,u,9){0)\\\l, (3.8) 

Jo 

where iyo> Q, C = C{eQ) are independent oft. 
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The proof if postponed to the end of Section 13.11 

Lemma 3.2. There exists 0 < eo < 1 and ho > 0 such that if E < e < eo and h < ho, then 
for a suitable /3 > 0, there holds 

||(p,M,s)(t)||^ + /3||Vp(t)f+ 1/0 [ II V(p, u,s)(r)f dr+ J /0 
<C\\\{p,u,0m\\\l + Ch'^ ds, 

Jo 

for some constant C > 0 independent oft. 

Proof. Under the transform of (Irel) . the system ( 12 . 11 ) is transformed into the following system 
for (p, u, s) 


||ftAp(r)fdr 


(3.9) 


dtp + u ■ Vp + (1 + p)divu = 0 , 


dtu^ 


V((l+p)^(l + s)) 


^ Am — ^ ^ VdivM + u ■ Vu 


dts + u ■ V9 — 


1+p P+1 P+1 

h^ AVp h^ div(VA/p + 1 ® Vi/p + 1) 

12 p -\- \ 3 p \ 

2n f Vs 2(l + s)Vp'| 

soTT) lT+7^ 3(i + p)2 / 

Vs 2(1 + s) 


—* 9h 


-K7(7- 1 ) 
2 

^3(1+ p)^ 

h? 


Vp\yp 


(l + p)2 3(1 + p3) 

{^|Vu+ (Vm)^P + A(divM) 2 | 

:div((l + p)Au) =: hn,. 


36(1+ p)T' 

Recall that E^{p, u,s) is given in (1X71) . We compute 

dtE^ =(1 + p)u ■ Ut + {— H-^~r(l + *((1 + pY ^ ~ 1)) + s 


+ 


7-1 

s^}pt + {^-r((i + pF - 1 - 7P) + 2-^-^}st 


2 ( 7 — 1 ) ' '■7—1 '' ' ' 7—1 

=div{- • • } — p|VmP — (p + A)|V • Mp — 2 h ;(7 — l)Vs • Vp 
- «:(7 - l)2|Vpp - AciVsp + 0{E)\D{p, u, s)|2 

+ (1 + P)u ■ gn. - {^-r((l + P)^ - 1 - 7P) + . 


- 7-1 


7 - 1 


(3.10a) 


(3.10b) 


(3.10c) 

(3.10d) 


(3.11) 


II 

Now, we consider the integration in space of the last two terms I and II on the RHS of 
(13.111) . For the first term /, by integration by parts, and using (I3.10al) . we obtain 

I Idx = -^Jv-uAp+^J Vm(Vp 0 Vp/(1 + p)). (3.12) 

The last term on the RHS is easy to be bounded by 
h? 


12 


Wu{Wp® Vp/(1 + p)) 


<c||VM||ioe||Vpr <cft^F;||Vp||' 
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For the first term on the RHS, we use (I2.1al) to obtain 


12 


V • = 


12 


12 


VV • itVp 
/■ 1 


1 + p 


V{pt+u- Vp)Vp 


d 

f |VpP 


f PtiVpP 

/ 

24 dt J 

1 + p 

24 j 

(1 + p)^ 

24 J 


1 + p 


'-Vp. 


But from (|2.1ap , it is easy to know that 

< C||div(pu)||L=o < CE^ < CE, 

and by integration by parts, 


V(u • Vp) 

1 + p 


Vp = - 


diU ■ Vp 

1 + p 

diU ■ Vp 


dip- 

dip + 


Ui ■ diVp 


1 + p 

V-u|Vp|2 


Vp 


J 1+p J 1+p 

<c'(||Vu|U= + Nloo||Vp|U~)||Vpf 

<CF;||Vpf. 

Therefore it is easy to see from (13.121) that 

|Vpp 


u -Vp 

( 1 +^ 


|Vp|' 


d 

Idx> — — — 
- 24 dt 


1 + p 


— CfrE\\D{p^ u, s) 


On the other hand, for the second term II in (13.111) . we have 


//dcc<5o||(Vp,Vs)|| 




IIAuf. 


In addition to (13.111) . we compute 


(l + p)2 1 „ (1 + P)"v. 2(1+ p) „ 

-u} ■ Vpt + 1 -——r^P • ■ ^ppt 


_a 

dt 


= Vp- 


2p + A 


=div{...}- 

(divu)^ 


1 


2p + A 


Vp • Vs — 


2p + A 


+ 0(£')|0(p,u,s)| 


2,11 “h A 

2 , (1 + p)" 


2p + A 


2p + A 


^p-gn- 


J 


After integration in space we obtain for the last term, 


/ 4 / '0b ^ + f ■ ■i™(v+TT® Vvi^rTT) 


— : Ji + J 2 . 

For the term Ji, we have 

Ji=- 


(1+^)|A |2 (l+f')|V7 |2 a 


/■ (1 + P) I 


> - 


127 2p + A 

^ f b + P) 
24 J 2 p + A 


2fi “h A 

|Ap |2 -Ch^Fl^llVpf. 


(3.13) 
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For the term J 2 , we have 


J 2 < 


f{i+p) 


+ CteE^\\Vp\\^. 


48 J 2/i -j- A 

Multiplying p.l3p with a constant /3 and integration in space, and then add the resultant 
to (13.111) integrated in space, we obtain 


d 

dt 


Wp ■ u] dx 


+ / p\'S/u\'^ + (/r + A)|V • + -kVs • Vp + -k|Vp|^ + Kl^sf dx 

J 3 9 


2p A 


1 


Vp • Vs 

— - 1 - 

2p -(-A 3 2p A 


|Vpp- 


(divu)^ 1 + P 


(3.14) 


+ 


2p A 48 2p A 

CK^. 


I Appda: 


<0(£;)p(p,«,s)|p + <5o||(Vp, Vs)f + . 

Oo 

Note that as in [19], if we take /3 small such that 

0 < /? < min|^^j^^-^,(p + A)(2p + A),4k(2p + A)| , 


where p 2 is given in Lemma l3.11 then 

l;°(p,«,s) + /3Q|Vp|2 + ^^ 




Vp-u) >l(p2 + Is2 + |«|^) + A|Vp|^ 


and 


p|Vup + (p + A)|V - Mp + hvs •Vp + h|Vpp + «;|Vsp 

3 9 


■/3 


Vp-Vs 5 1 1^ |2 (divii)^^ ^ ,9 . 5k, 


2p, X 3 2p -|- A 


|Vp|^- 


2/i “h A 


>p|Vup + —|Vsp + -«;|Vpp. 


Integrating in time over [0,t] and taking Sq and E sufficiently small (say, Jq = 1/20), we 
obtain 

f i + ^|Vsp + |k|VpP + J App 

^ Jo ./R3 I 


27' 


3 

<[ ■E'’(p,m,s)+ /3 (i|Vpp + Fl^Vp-A da; + CK^ f ||Au|pds 
JR3 V ^ + A / Jo 

<C|||(p,u,P)( 0)|||2 + C;.4 /‘llAzrfds. 

Jo 

Now, properly choose the constant t'o and C > 0, we finish the proof. 

Lemma 3.3. Under the same condition in Lemma \3.‘A we have 

{||Vu(t)|p + ?i^||Ap(t)|p} + f ||IiAM(s),fiVV • u(s)||^ds 

Jo 

<Ch^\\Vuo,hApof+ Ch^E [ \\Vpfds + 6if || VP, Vpfds + ^ / ||Aufds, 

Jo Jo 0l Jo 

for any positive constants (5i > 0 and t > 0. 


96 2p + A' 


(3.15) 


□ 


(3.16) 
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Proof. We now take the inner product of (I2.1bl) with —hfAu to obtain 


A d f, „ „ 


I Ai 


1 + p 


+ {^l + X)h^ 


(VV • u) 

1 + p 


■ Au 


=hf j{u ■ \7u)Au + y y 


' + 1 


f div(VVp + 1 (8) Vy/p + 1) 

T J rTp 


P+1 

5 

Au =: E] di. 


Vp • Au — — 


f AVp 


12 7 1 + P 


Au (3.17) 


For + 3 , we have by integration by parts twice 
,2 f |VV-U |2 


+3 —(p+ x)h^ y 

- (p + \)h‘^ 
>(p + A)p2 / 


(1+p) 
f dip 


I ( 1 +p)^ 

VV-u'2 


i9fcV • ui9feu® + (p + A)P^ 


9fcP 


(1 + p) 


rSiV • udku^ 


(1+p) 


-CP^F;||Vp, Au||^ 


since ||Vu||l“ < ClluHija < CE. 

For i?i, after integration by parts twice, we obtain 

Ri = -ti^ [(u • Vu)Au - [ dku^dku^d^u^ + / V • u|Vu|^ 


which implies that 




Ri=-ti^ J dku^dku^diu^ + y y ^ ■ 

For i ?2 and i? 3 , by Holder inequality we obtain 


(jfp 

R2 + R3< di(||VP, Vpf + —IIAuf, Vdi > 0. 

Pi 


For i? 4 , we have by integration by parts and (I2.1a|) 


i?4 — — —— 


h* f ApVp 


12 J (l + p)2 
r ApVp 
uj (l+p)2 


• Au + — 


f Ap 


12 J 1 + p 


AV • u 


• Au 




Ap 


j-{9tAp + [A, 1 + pjdivu + [A, u]Vp + u • VAp} = E din- 


12 J (l + p)2 

It is easy to show the following estimates 

i?4i < CP^IlVpIlioollAp, Auf < C'F;ft2||PAp,PAuf, 


hf d 

ddi2 = ~ 777 TT 


|App 




dtP 


^ d 


2AdtJ (l + p)2 12 J (l + p)3 

|App 


|App 


24 dt J (1 + pY 


+ CEhY\fiAp\\Y 
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thanks to (13.31) and 

i?43 + Raa <C'?i'‘||Ap||(||[A, 1 + p]divM|| + ||[A,u]Vp||) 
<C;i4||Ap||(||Au||||Vp|U=. + ||Ap||||Vii|U.) 
<CEh^h/\p, fiAwf, 
thanks to Lemma 12.61 and by integration by parts 


= 


12 J (l + p)2 

f |Ap |2 




-2AJ (1 + p) 2^^^“-T 


12 7 (l + p)2 

|App 


div.-^^ 


|App 


6 7 (1 + p)^ 


u -Vp 


i^ + pY 


u - V p 


<CEh^\\hi^pf. 
Therefore, we obtain 


hY d 

Ra<- 777-r: 


|App 


+ CT;;i2||ftApf. 


24 dt J (1 + pY 
For the term R^, it is easy to show that 

i?5 <CEh^\\hApf. 

Hence, putting all the estimates together, we have from (I3.17|) that 

|App , ,2 f IA11I2 


f 2 _ 

2 dt J 24 dt 7 (1 + p)^ 
<CtYE\\yp,Auf + 6Y\y{p,0)f 


ph 

CtY 


1+ P 


(p + x)h^ 


|VV-w|' 


(1 + P) 


II Aitf + CFlti^lltiAp, tiAuf. 


Take £0 and Hq small, then for any e < Eq and h < Hq, integration in time over [0,t] yields 
the result for any positive constant (5i > 0. □ 

Proof of Proposition of \S.l[ By (13.6|) and Lemma IXTl the left hand side of (13.91) is equivalent 
to the norm 

II(p, u, e){t)f + 711 Vp(t)f + uo f II (Vp, Vw, V0, hAp)is)fds, 

Jo 

for some another suitable constant uq > 0. Hence (j3.9p implies that 


||(p,w,0)(t)f+7||Vp(t)f+ U0 / ||(Vp,V«,VP,tiAp)(s)fds 


<C|||(p,«,0)(O)|||^+C';i" / IIAixll^ds. 


(3.18) 


Now, taking and Hq in p.l6l) sufficiently small, say, = uo/4 and < vq/ACeo, we then 
obtain 


;i2{||vu(t)f+ ;i2||Ap(t)f}+ / \\{HAu,hyy■u){s)fds 


<CfY\\Vuo,fiApo\\ 


2 , ^0 


||(Vp,V0)fds- 


4Ch‘ 


4 /-t 


(3.19) 


Jo 


II Au|pds. 

















12 


X. PU & B. GUO 


Add (|3.18l) and (I3.19|) together, and then taking Hq even smaller such that C/ig + dC/iQ/uo < 
1/2, we then obtain 

|||(p,n, 0 )(t )|||2 + i.o j\\{Vp,Vu,Ve,hl\p,hl\u){sWds <C\\\{p,u,e){ml (3-20) 

Jo 

for some positive constant ug > 0 depends only on p, A and k. in particular, ug and C are 
both independent of t. □ 


3.2. Higher order estimates. In the following, we denote 5“ = the partial 

differential derivative operator with multi-index a = (01,02,03). For our purpose, |o| < 3 
suffices. We sometimes abuse the notation to use o ± 1 to stand for o ± /3 for a multi-index 
with |/3| = 1 and o > /3 in the case of o — /3. We will prove the following 


Proposition 3.2. Let a be any multi-index with 1 < |o| < 3 and s = |o|. There exist some 
constants erg > 0 and Hq > 0 such that if E < Eq and h < Hq, then the following a priori 
estimates hold for all t € [0,r], 


|||d“(p,u,6»)(t)||p-h ug / ||5“V(p,'u,6», Wp, WM)(r)f dr 
Jo 

<C|||5“(p,u,6»)(0)||p-hC f \\V{p,u,0)fH,-idT + C f \\hA{p,u){T)fHs-idT, 

Jo Jo 


(3.21) 


for some ug > 0 and C = C'(eg) independent oft. 


This proposition is proved as a direct sequence of the following lemmas. 

Lemma 3.4. Under the assumptions in Provosition 1 3. SI there exists some constants dg < 1 
and Eg < 1 sufficiently small, such that 

\\d°e{t)f + K[ ||d“vd(r)fdT 

Jo 

<||5“6>o|p-h(5g(p-hA) [ ||d“V • u(r)|pdr-h [ ||d“AM(T)|pdT (3.22) 
Jo ^ 0 ^ Jo 

+ Y f l|V(p,n,6>)(r)||^,-,dr-H ^ / \\Au{T)\\^^^_,dT, 

do Jo doK Jo 

for all d < dg and E < e < eg, where C is independent oft. 


Proof. Applying d°‘ to (I2.1c|) and then taking inner product of the resultant with d°‘6 to 
obtain 


l_d 
2 dt 


Wd’^er 


2k 

T 


|d“V6l|2 

(1 + p) 


2k 

T. 


Vd“6i • Vp 

(1 + p)^ 


a“6» 


- J a“(w • V6i)d“6» - ^ y ((6» 4- i)v ■ “) + ^ / ' 5 “ ( 


/div((l -f p)Au) 


0 + p) 


d“6» 




fj,\Vu + (V'u)^p + 2A(divw) 

(i+p) 




For the first term i?i, we have 


Ri <CK\\d°‘ve\\\\\/p\\L3\\d°‘e\\Le < c«r£;||d“V6»f. 
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For the term we have 

i?2 <||5 “(m-V0)|| 

<(||w|U3||a“V0|U2 + ||9%|Ue||V0||i3/2)||V5“0|U2 
00 K 

thanks to Lemma 12.61 For the term i?3, since |0 + l|<3/2by (13.31) , we have 


i?3 <(||0 + 1 ||l-||9“V • u\\l^ + ||a“0|Ua||v • u\\L^)\\d^e\\L^ 

<5oK\\vd^e\\l, + 5^{^l + A)||a“v • u\\l. + 

Oo 

For the term i?4, by integration by parts, 

<Ch^d<^+^9\\ |||9“-Miv((l + p)Aii)|| + ||[9“-i, Y^]div((l + p)Aii)||| 

<C'h2||a“+i6»|| {||a“-MivAM|| + ||Au||L~||5“-Vp||i2} +C'h2||a“+i6i||x 

{l|V(^)||L<»||a“-2div((l + p)Aii)|U2 + ||^||^._,,e||div((l + p)Aii)|U3|, 

where s = |q;|. Then making use of Lemma 12.51 and 12.61 and (I3.1D - (I3.5L one obtains 


i?4 <Ch^d-+^e\\ (IIAull^. + IIAuIU^IIpII^.) 

+ Ch2||a“+id||||Vp|U..(||Au||^,_, + ||p||^._,,e||Auy^ 
+ c;ila“+id||||p||^._,.a(||vAu|U3 + \\vp\\L<^\\Auh^) 


<do«:||9' 


ia+1/ 


+ 


Ch'^ 

Sqk 


IIAull 


CE^ti^ 

Son 


OqK 


Similar to i?4, we have for the term that 


D 2 /■ /^/i|Vu + (Vu)^|2+2A(divu)2 

''== 3 /^ i- JTp) - 

<3„K||a“+‘ef + ^||(vp, 

OqU 

Putting these estimates together, we obtain 


+ yJ + ^o(m + A)||a“V• uWh 


c 

<5o 


wor 


L2 


Cff 

5ok 


WAuf + ^||(Vp, + ^^||Auf^._,. 

0{)rx 0[)l% 


Integrating this inequality in time over [0,t] and noting 1/2 < 1 + p < 3/2, we know that 
there exists some constants i5o < 1 and eo < 1 sufficiently small, such that (j3.22|) holds for 
all d < Ao and E < e < eq. □ 
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Lemma 3.5. Under the assumptions in Provosition lS.Sl there exists some constant £o < 1 
sufficiently small and Hq < 1, such that 


\\d°^[p,u,hS/p){t)f+ vq / \\d°^Vu{T)fdT 


<CW{p,uffiVpmf + CE f \\V{p,u){T)\\lffiT + C f \\Vffiuffi){T)\\\^_ffiT. 

Jo Jo 

for all E < e < Eq and h < fig, where C is independent oft. 


(3.23) 


Proof. Applying 9“ to (I2.1bp and then taking inner product of the resultant with 9“tt, we 
obtain 


1 d 


Au 




, V divM 


)d° 


P+1 J P+1 

= - j d°‘{u-Vu)d‘^u- J d°‘VBd‘^u- J 9 “(^^Vp)( 9 “u 

, /’ oa ^^ VP^aa /■ s , a / div {---} 

+ 12/^ ^ = 

Now, for the term L 2 , we have by integration by parts 

f d°‘Au f 1 

+2 = — P / - d°‘u —p / [i9“,-]Au9“u 

Jl+P J 1+P 


=P 


|a“Vu|2 f 9“Vu 




1 + P "”7 (1 + p)^ 

Invoking Lemma 12.51 and 12.61 we obtain 

l|[9“,-^]Au|La/B <C||Vp|Ua||9“-iAu|U2+C||Au|Ua||5“( 


1 + P 


1 + P 




<C7L;||a“-^Au|U2 + CE\\Vp\\Hs-^. 


(3.24) 


Hence 

+2 >^||9“Vuf - Cp||a“Vu|U2||Vp|Ua||a“u|Ue - C||[5“, ^]AuLa/5|ia“u|Ua 
3 1 + P 

> YP“Vuf - CE{1 + p)||5“Vuf - Cii;||Vp||^._, 

>|||5“Vuf-CA||Vp||^._, 

by taking A < p/6C{l + p). Similarly, we have for L3 that 

+3>^lia“V-uf-CA||Vp||^._, 

For the RHS term Ai, we have by integration by parts that 

Ri =i y V • u\d°‘u\^ - J[d°‘,u]Vud°‘u 

^CWVuUYld^^uf <CE\\Vu\\l^_,. 

For the term R 2 , we have for any dg > 0 that 

R2<Sop\\vd^uf + -^\\ve\\%_,. 
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The term R 3 will be treated with much more effort, from which some good terms will appear. 
By integration by parts, 


R.=- 


p+1 J p+1 


= f ^ii-5“pa“divu+ f f = Vi?3z. 

jp+i jp+i J p+i ^ 

It is easy to show that 


i?32 < C||V(0,p)|U~||a“(p,u)f < CE\\Vip,u) 


Ijys-l, 


and 


i?33 <c||a%|| (||v(^)iu~||vp||^._, + ||vp|u^||^ + 


i + p 




I + P 

<CE\\vip,u,e)\\%_,, 

thanks to Lemma I^TSl Differentiating the continuity equation (12.lap with 9“ yields 
(1 + p)9“divu = -dtd‘^p - a“(u • Vp) - [5“, 1 + p]divu, 
which implies that 

i? 3 i = - y j^^^d°‘p{dtd°‘p + 9“(u • Vp) + [5“, 1 + p]divu} = ^ Rsu. 
It is immediately from (j3.3p and (13.41) with p = 3/2 that 




dt 


' + 1 


■Wp? 


(i + p)= 




9 + \ 




2 dt J (1 + p)^ 

0 + 1 


\dy\^ + c\\dte,dtp\\L^f4dy\\h 

|a“pp + C'£;||9“Vpf. 


2 dt J (1 + p)^ 

For the term R 312 , we have by integration by parts that 

<q|v(p,u,0)|k<.||a>f+ ||a>||(||Vu|U<.||Vp||^._.||Vp|U~H 

<CF;||v(p,u)f^._,. 

The same estimate hold for Rsis- Combining all the estimates for R 3 , we obtain 


i?3 < — 


1 d 


9 + 1 


)9“p|2 + cF;||5“Vpf + cf;|| V(p, u, ( 


2dtJ (l+p)2'^ '-'^'1 - v^„ , 

Now, we consider the estimate of R 4 . By integration by parts, we obtain 




diP 


Ra =j^ I - TV / ^“ ( 


> +1 


12 


1 


P +1 


]9,p 


4 /a‘(;;?T)3“Vdiv„ + ^/3- ([Vt^iv,) 

=— / -^9“Vdivu+— [d°‘, -]Vpa“Vdivu 

12 7P+1 127 p+1 


+ 


12 


[V2, 


p +1 


]Vp)a“+iu=:^i?4.. 
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Using the continuity equation (I2.1al) and similar to the term R31, it can be shown that 
ii’41 = - -^ / {dtd^^^P + 5“V(u • Vp) + [a“V, 1 + p]divM} = ^ i?4H. 


For i? 4 ii, we obtain 


< - 

By integration by parts, 

4 /<!- 

24 ,/ 


d r 

|5“Vp|2 

24 dt J 

(1 + 4 

d r 

|5“Vp|2 

24 dt J 

(1 + 4 

9“Vp 

( 1 + 4 '' 

• a“v 4 


-Cn^E\\d°‘Wpf. 


12 J (l + p )2 


0 - u / 


(i + p) 


[,9“V,w]Vp 
[9“V,u]Vp, 


,( 1 +P)^ 

and hence by Lemma 12.61 and (13.21) 

R 412 <cn'^\\v{p,u)\\L’^\\d‘^^p\\^ 

+ c;i2||9“Vp||(||9“Vp||||Vu|U~ + ||Vp|U^||a“VM||) 
<Ch^E\\d’^V{p,u)f. 

Similarly, by Lemma HH] and (13.21) . 

i?4i3 <creE\\d<^vip,u)f. 

For the term R 42 , we have 


i ?42 = - ^ [ v([9“, ^-]Vp)a“divM 

127 p +1 


12 


/^]V49“divu 


ip+ir 

But by the commutator estimates, we have 

7-^]VpIIl^ <l|Vp|4.-2,a|l7-^llL3 + IIVpIU^II- 


4 + 4 


4 + 4 


4 + 1 ) 


2 II 


<cf;||Vp||^., 


and 


hence 


<C£;||Vp||^., 


Similarly, for R43, we obtain 


i?42<Cn4||V(p,u)||^.. 


i?43<C;i4||V(p,u)||^.. 


Putting all the estimates for R 4 together, we obtain 

d f |a“Vp |2 , 


i?4 < - TTT^ 


2AdtJ {l + pf 


Ch^E\\V{p,u)\\ 
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Finally, for R^, it is easy to show 


/■;,a-i/div{Vp(g)Vp/(l + p)}^^„+i 


^5 =Y^ / 9“-'( 


P+1 


-)d° 


<Ch^E\\d‘^yip,u)f + CfeE\\yp\\l^_,. 

Now, putting all these estimates for (13.241) together, and taking Jq = 1/4, we obtain, 

1 d 


2 dt 


119' 


c.,„2 , 1 d f 0+1 d /-IS^Vpl^ 


2 dt J (1 + p)2 


2AdtJ (l + p)2 


+ ^||9“Vixf+ ^||9“V.nf 


<C\\d^9\\^ + CE\\d‘^Wp\\^ + CE\\V{p, u, 9)\\%^_, + Ch^EWd^^Vur 
Integrating in time over [0,t] completes the proof, thanks to h < 1. 


□ 


Lemma 3.6. Under the assumptions in Provosition 1 3 . ill there exists some constant Eq < 1 
sufficiently small and Hq < 1, such that 

h^\\Vd°‘(u,pffiVp){t)\\'^+ VQh? f ||A9“M(T)|pdT 

Jo 

<Ch‘^\\Vd°‘iu,p,h'Vp){0)f+ Ch^E [ \\Ad°‘u{T)fdT + Ch^E [ ||A9“p(r)fdr (3.25) 

Jo Jo 

+ Ch^ [ Wd^^Wdiffillffir + Ch^E f \\V{p,u,9){T)\\lffiT. 

Jo Jo 

for all E < e < Eq and h < fig, where C is independent oft. 

Proof. Applying 9“ to (I2.1bl) and then taking inner product of the resultant with —ft^A9“M, 
we obtain 

+ d^{^)Ad^u + {p + X)ffi 1 9“(^^)A9“« 

J 9“(w Vw)A9“u + ft2y’ 9“V6»A9“u + y p)Ad°‘u 




- ^ / 9“(^)A9“ix + ^ i?, 


div{...}, 


127 P+1 3 

Now, for the term + 2 , we have by integration by parts 

qa,,|2 


P+1 


+2 =ph^ J + Pli" J [5“, Y^]AuAd^u. 


Since 


1 


||[9“, ^]Am|U= <C||Vp|Uoo||9“-iAu|U. + C||Air|U3||9“(-—)|Ua 
1 + P 1 + P 


<CE\\Au\\^..,+C\\Au\\h.\\p\\^.,, 

<CE\\{p,u)\\^,+^, 

thanks to Lemma 12.51 and (1^ . we have 

L2>^\\Ad‘^ur-CffiEmp,u)\\%^^,. 


Similarly, for the Lg, we have 

(p + X)ffi 

J^3 A -- 


II9“VdivufII (p, u) 


lU- 


2 
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For the RHS term i?i, we have by integration by parts twice that 

Ri = -h^ J djd°^{u^ ■ d^u)djd°‘u 

= — J u^djd°‘diu ■ djd°‘u ~ J [djd^^^u^jdiudjd^^u 

=y J d,u^\djd°‘u\^ - J[djd‘^,u^]d^udjd‘^u 

<Ch^E\\u\\l,^,. 


For the term R 2 , we have 

R2 = h^ / d°‘VeAd°‘u < donh^\\Ad°‘uf + —||a“V6»f. 

J OofJ. 

By integration by parts, 

R-i=-h^ [ V5“(^-^Vp)Va“M 
J P+1 

= -h^ f ^Jllv9“VpV5“u - f V9“m[V9“, ^-^iVp 
J P+1 J P+1 

=n'^ f i±lva“pVa“divM + f v(^-^)V5“pV9“m 
7p+1 Jp+1 

-h^ /'v5“M[V9“,^]Vp = y i?3.. 

J p+1 ^ 

It is easy to show that 

R 32 < h"||V(0,p)|Uoo||a“V(p,n)f < CteE\\Vd^{p,u)f 


and by Lemma [2+] and [2Rl 

R 33 <Cft2||V5“n|l||[V5“, ^]Vp|| 

P+1 

<Ch2||V9“u||(||Vp||^.||V(p,0)|U=o + ||Vp|U~||(p,0)||^.+O 

<Ch^E\\V{p,u,0)\\%. 

Differentiating the continuity equation (I2.1al) with 9“ and then inserting the resultant to 
i? 3 i, we obtain 

R 31 = - ^" / ■ ^P) + 1 + ^ R3U- 

It is immediately that 


R 311 = — 


d 

2 dt 


^ d 
~ 2 dt 


8 + 1 
(1 + p)^ 
8 + 1 


|V5“p|^ 


2 


5*(7r^)|V9“pp 


(i + p)^ 

\yd°‘p\^ + Ch'^E\\d°‘Ap\\^, 


( 1 +p)^' 

thanks to (1531) and (13.41) again. For the term R 312 , we have by integration by parts that 


«3.4/div(('’ 


l)u 


( 1 +p)^ 
<Ch^E\\Vd‘^{p,u)f 


)\yd^p\^-h^ I ^^^V5“p[V9“,n]Vp 
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The same estimate hold for i?3i3. Now, we consider the estimate of R 4 . By integration by 
parts, we obtain 


h* f Ad°^pAd°‘u_ 


-Aa“divu-,„ 
12 y p +1 12, 


(p+l)2 12 




= ^ i?4i • 

For the last two terms, it can be shown that 

R42 + R43 <Ch'^\\Vp\U\Ad‘^{p,u)f + Ch'^\\Ad'^u\\\\AVp\\L4d‘^{l/{l + p))\\lo 
<Ch'^E\\Ad^{p,u)\\^+ Ch^E\\yp\\l,. 

Using the continuity equation p.la|) and similar to the term R 31 , it can be shown that 
h'^ f Ad^p 


-R41 — — — 


< - 


12 7 (1 + p)^ 

d f |5“Ap|2 


2 AdtJ (l + p)2 
Finally, for i?5, we have 


{dtd^Ap + d°^A{u ■ Vp) + [a“ A, 1 + p]divu} 
Ch^E\\d‘^A{p,u)f. 


R, / a"(AU^^^Ai±dl)Aa". 

127 p+1 

<Ch'^E\\d‘^ A{p,u)f + Ch^E\\v p\\\^. 

Now, putting all these estimates together, and taking 5q = 1/4, we obtain, 


y|l|va“uf+ u^l|Aa“iir + 


a„,l |2 , , (^ + ^)^^l oa II 2 


||a“Vdivu||^ 

2 


+ 54/A±i,|v3VP + ^d/'| 8“^'-l 


2 dt J (1 + p)2 


2AdtJ (l+p)2 


<crF;|| A9 “m||^ + ch^E\\Ad°‘pY + Ch^WveY + C'ru||v(p,u,6»)||^,. 

Integrating in time over [0,t] completes the proof. 


□ 


Lemma 3.7. Under the assumptions in Provosition I S. SI there exists some constant £o < 1 
sujficiently small and Hq < 1 , such that 

(2p + X)/3\\d-^Ap{t)f + /3 r ||5'^A(p, hVp){T)fdT 

Jo 

<C'/3||5''Apof + /3||a^VMof + Cl3\\d'^divu{t)f 

+ C/3 f Wd'^AOf dr+ Cfd f Wd^^VuWl.dT + Cp f \\V{p,u,9)\\l^_,dT, 

Jo Jo Jo 

for any P > 0, E < £ < Eq and H < Hq, where C depends only on p, A and some Sobolev 
constants. In particular, C does not depend on h > 0 or t > 0. 
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Proof. Let 7 be a multi-index such that I7I = |a| — 1 and 7 < a. We apply d'^ to (I2.1b|) and 
then take the inner product of the resultant with Ap to obtain 




-13 J d^ut ■ 5'^VAp + I3pjd'^ (1^) 

+ i3{p + X) J d'^ (tT^) J 

- /3 y d^VOd^VAp J d^VAp 




Q1 


/ VAp 
U + P 


d'^S/Ap - 


P?P 


97 


( dW{---} 
V 1 + P 


d^VAp = 0. 


(3.26) 


We first note that by Holder inequality 


L 7 


PfP r |5^VAp|2 

~~[2 J 1 + p 

>^lia7vApf- 
>^||5'^VApf - 


c/3?i2||a7vAp||(||vAp||^._.||Vp||i^ + ||vAp|U3||a'^(^)|Ua) 

CptPE\\Vp\\l^-CpfpE\\Vp\\l^_,. 


and by integration by parts and Holder inequality 


Lq =P [ d'^Apd'^Ap + P [[d'^'V-,—-—Wpd'^Ap 
J 3 + p J 1 + p 

>^WApr-CPWAp\\ (^||Vp||^._.||V(p,P)|U=o + ||Vp|U~||5^V(^)||) 
>-^\\d^Apr - cpE^yip,9)\\l,_,. 


For the term Lg, we have 


\Ls\<SiPfP\\d'^WAp\\ 


2 , Cph^E\^^ „2 




-livpil 




For the term L5, we have by integration by parts 


\L.\ = 


P J d^AOd^Ap 


</3<5i||a'^Apf + J-||a^A0f. 


For the term L4, we have by integration by parts 


IL4I = 


P 


Vd'^{u-Vu)d'^Ap 


<P\\d^Ap\\{\\u\\L^\\yu\\Hs + ||Vn||^._,||Vw||L~) 
</35i||9^Apf+ ^||Vu||^._,. 
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For the term L3, we have by integration by parts, 


(m + A)-1L3 = - /? / • (^^)d^Ap 

J 1 + P 

= - 13 j Ap - 13 Aa''V, Y^i vdivtia''Ap 

^ a^A((i + rtdiv..) g, ^ r [a'-A.i + dldiv~ g, 
J (1 + P) J (1 + P) 

- ;9 y [9'^V, Yy-^]Vdivu5'>'Ap =: ^31 + L32 + ^33- 
For the first term ^31, using (I2.1al) . it is easy to show by integration by parts, 


/3 d /3 

2dtJ + (l + p)3 2 


+ 



~ 2 dt 


/ 


[a'>'A,u]Vp 

(1+p)^ 

I^T-App 

(1+p)^ 


d^Ap 

C/3i3||a^A(p,ii)|| 


2 


div( 


U 

(1+^ 


)|5'^App 


thanks to (13.31) with p = 00 and Lemma 12.61 By commutator estimates, 


|3>32| <C'/3||a3'AM||(||Vp||Loo||divu||^, + ||a3'Ap||||divM||L=o) 

<CI3E\\^iP,n)\\%, 

|i33| <C/3||a'^Ap||(||Vdivn||^..,||Vp|U.» + ||Vdivw|U3||p||^.,e) 

<C( 3 E\\^iP,u)\\l.- 


The term L 2 can be treated similarly and will lead to 

Pp d f \d'^App 


L 2 >- 


-C/3F;||V(p,n)rp,. 


2 dt J (1 + p)2 

Now, we focus on Li. Integrating Li in time over [0,f], we have 

/3 f Li{s)ds =P f f d'^davutd^Apdr 
Jq Jo j 


/o 

= P [ d'^divud'^Apdx 
Jr3 

= : I(t) - /(O) - /3 f Li2{s)ds. 
Jo 


0 Js.^ 


diYud'^ Aptdxdr 


By direct estimates, we have 

\I{t)\ <di^||a^Ap(f)||2 + ^\\d^divu{t)\\^, 

Oi 

\im<PWApof+P\\d^yuof. 
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For /3 fg Li 2 (s)ds, we have by integration by parts and (I2.1a|) 

—/3 f Li2{s)ds = — j3 [ f V9’*'divuVi9^(it • V/9)dxds 
Jo Jo dR 3 




0 


Vi9'’'divuVi9'’'[(l + p)divu]dxds 


<CP / ||V5^divu||(||u|U^||Vp||^. + ||^r||^.,a||Vp|U3)dr 


+ C/3 / ||V5'>'divu||(||divu||^. + ||divu||i 3 ||p||^,, 6 )dr 
Jo 

<CI3 [ WVd^dwufdr + CpE^ [ \\V{p,u)\\%Jt. 

Jo Jo 

Now, we fix some (5i small, depending on p and A, such that 5i = min{(2^ + A)/18,1/142}, 
where C is the constant appearing in the above estimates. Then for such (5i, fix some small 
£0 such that £o < min{l, 1/(72C')}. By integrating ()3.26p in time over [0,t], and then taking 
E < So small, we complete the proof. □ 

Now, we prove Proposition [221 

Proof of Provosition [STR Adding the estimates in Lemmas 13.4113.71 together. First, we take 
<^o = uo/4(^ + A), and then hg small such that hg = min{i5oUoK/4C, k/AC, 1}, then choose /? 
small such that CfJ < CfJg := min{K/4, i/o/4,1/4}, and then for such fixed /3, we choose eo 
small such that Ceg < min{z/o/4, k/ 4,/3/4}, then there holds, 

W{p, u, 9, Wp)(t)f + h^Wd^iu, ftVp)(t)f + /3||a^Ap(t)f 

+ uo / |15 “V(m, 0, W'u)(r)|pdr+ /3 f ||5^A(p, W/9)(r)|pdT 

Jo Jo 

<c||a“(p, u, 9, ftVp)(o)f + va“(u, ftVp)(o)f + c\\d^Apgf 

+ C [ \\Vip,u,9)\\l^_,+Ch^ [ ||V(p,u)(r)||^,dr, 

Jo Jo 

for all ^ < ^0 7 where C is independent of t. Rephrasing this in the ||| ■ ||| norm, we obtain 


|||5“(p,u,0)(t)||p+uo / ||9“V(p,'u,0, Wp,ftVu)(T)f dr 

<^1 P“(p, u, 0)(O)|| p + C f II V(p, u, 0)||^._, + Ch^ f II V(p, u)(T)||^.dr. 

Jo Jo 

The proof is complete. 

Now, we are ready to obtain a priori estimate for the solution of m- 


□ 


Theorem 3.8. Suppose that for some T > 0, {p,u,9) G £ 3 ( 0 ,T) is a solution of (12.1|) 
satisfying E < maxo<t<T lll(0)'w,0)('r)|||3 < e. Then there exists some eg > 0, uq = vg{eg) > 
0 and Cg = C'o(£o;Uo) such that 

I||(p,u,0)(t)|||2 + Z2o / ||iA(p,u, 0 ,ftVp,ftVu)(r)|| 2 , 3 dr<C|||(p,u, 0 )(O)|||i 
for any {p^u^9){t) satisfying E < Sq. 
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Proof. Add (13.211) for |q;| = 1 to (^ + 1) times (13.81) to obtain 

\\\{p,u,9){t)\\\l + ,yo [ P(p,«,0,ftV«,Wp)(s)||Uds<C'|||(p,u,0)(O)|||i (3.27) 

Jo 

Then, we add (13.211) for |a| = 2 to (^ + 1) times (13.271) to obtain 

\\\{p,u,9){t)\\\l + ,yo [ p(p,«, 0 ,ftV«,Wp)(s)||^ 2 ds<C|||(p,u, 0 )(O)|||^ (3.28) 

Jo 

Again, we add (13.211) for |a| = 3 to (^ + 1) times (13.281) . to obtain 

\\\{p,u,9){t)\\\l + ,yo [ p(p,«,0,ftV«,Wp)(s)||^3ds<C'|||(p,u,0)(O)|||i, (3.29) 

Jo 

which is the desired estimate, completing the proof. □ 

Now, we prove Theorem 12.21 

Proof of Theorem \2.S[ The proof is easy by combining the estimates in Theorem 13.81 conti¬ 
nuity method and the local existence result. □ 


4. Proof of Theorem 12.31 

4.1. An improved estimate of Matsumura and Nishida |19| . The following estimate 
was obtained in m- Under the smallness assumption. 


\\{p,u,9)ml + ’^o / \\Dpir)\\l + \\D{u,9){r)\\ldr < CoH(p, u, 0)(O)||i 
^0 


and the following density estimates 

\\D*p{t)r - C\\DMt)r + ^0 f PV(r)f dr 

^0 


<C\\p4l + C\\u,\\l + C / \\DHn,e){r)r + \\D{p,u,9){T)\\ldT. 


(4.1) 


(4.2) 


Adding (IT^ to (C - 


G_ 

VO 

i4„^3M|2 


1) times (14.11) . one obtains 


\\{p,um)\\l + Wp{t)r + i^o / \\D{p,uMT)\\idr<c4{p.u,e)ml + c\\p4i ( 4 . 3 ) 


4.2. Proof of Theorem 12.31 To clearly specify the dependence of the solution on the 
parameter h, we denote {p^,u^,9^) the solution of the system (12.11) and (p°,u°,0°) the 
solution to (12.31) . First of all, from Theorem 12.21 and (14.31) . we have the estimates 


\\{p\u494mi + pV(t)f + h^ip^u^mwi + 

<C'(ll(Po,'«o,do)( 0)||3 -I- ||(po,Mo)(0)|l4 + IIpoIIs)- 


(4.4) 


and 


||(p^«^0“)(f)||^ + pV(t)r <CII(po,«o,do)ll^ + q|polli). 


(4.5) 


Now, we let 


N = p^-p°, U = u^-u°, e = 9^- 9°. 
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Then {N, U, 0) satisfy 

dtN + V((l + p°)U + Nu°) = 0, 


dtU- 


p'* + 1 


-AU- 


+ 1 + 1 


= • VC/ - [/ • - V0 - 


0 


Au° - 4^VdivC/- 
p'^T 1 

Vp^ - - 


p^“ + 1 ^ p° + 1 

AVp^ h"^ div(V a/p^ + 1 0 Va/p^ + 1) 


/ p + A 

P + A ^ 

\p^ + l 

po + iy 

^6»o + l 

e° + i\ 

^p'^ + l 

p° + 1/ 


Vdivu° 

Vp° 


9t0- 


12 pft + 1 3 

2k 2k 

3(TTp^^" ~ T 


1 


p'* + 1 
1 


1 + p^ 


-^{9^ + 1)\7 ■ u - ■ u 

o o 


1 + p° 


A9° = -w° • V0 - C/ • V6»° 
div((l +p'')Au'^) 


3(1 


2 

+ 3 


p" 

1 


36(1+ p'^) 

{^(V(u'' + u°) + (V(u'' + M°))^)(VC/ + (VUf) + A(div(M^ + M°)divC/} 
^ ^ { 2 (Vu°)^P + A(divu°)2} . 


1 + p^ 1 + p° 

Now, we multiply (14.6|) with N, U and 0, respectively, integrate the resultant over and 
then sum them up to obtain an energy inequality. Among the many terms, we only treat 
the following three typical thems in the following. First, for the viscosity term, we have 


>f J |VC/p-C||C/f, 


VC/C/ 


where the constant C depends on p and the norm of p^. Secondly, for the last term on 
the left of (I4.6bl) . we have 


p + A p + A 

pfi _l_ 1 pO + 1 


Vdivu°C/ = 


(p + X)N 


(p'* + l)(p0 + l) 
<C||Vdiv«0|Ua||A|U.||C/|Ua 
<g||VC/f+ C||(A,C/)||2„ 


Vdivu°C/ 


where the constant C depends on p and the norm of (p^, p^,u^). Thirdly, for the second 
to the last term on the RHS of (|4.6bl) . we have 


12 


< ^\\ur+^wp^rn. < ^iic/f+ ch\ 


where the constant C may depend on p and the norm of p^. The other terms, either 
depending linearly on the difference A, C/ or 0, or depending on the small parameter , can 
be estimated similarly. Therefore, we finally obtain after long but standard estimates 

~\\{N, u, emwl. + u|| V(C/, 0)f < C||(A, C/, 0)(t)f + Ch\ (4.7) 

where u depends on the parameters p and k, and C depends on the norm of || (p^, u^, 9^)\\. 
Similarly, taking inner product with A{N, C/, 0), one can obtain 

i|||V(A,C/,0)(/)||i. +u||A(C/,0)f < C||V(A,C/,0)(t)f+ 


(4.6a) 


(4.6b) 


(4.6c) 


(4.8) 
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where v also depends on the parameters and k, and C depends on the norm of 

\\{p^,u^,6^)\\. In the derivation of this inequality, we have the following term, 

/ 36(1+ + p'‘)Azi'‘)A0 < J|A0f + C;i^||div((l + p^)Au^)r 

<j\\Aer + ch\ 

where C depends on k and the iJ^-norm of p and 0. This procedure can not be proceeded 
into the higher norms, since the estimates then depends on the i/^-norm of u^. 

Combining the two inequalities (14.71) and (14.81) . we obtain 

i^lKiv, [/, 0 )(t)|| 2 ,, + H|v(t/,0)11?,, < cilKiv, u, 0 )(t)||?,, + C2h\ 

which implies, thanks to the Gronwall inequality, that 

||(7V,t/,0)(f)||?,, < [c2e^^Vci] 

In particular, we note that c\ and C 2 are independent of h. 

We also remark that we can improve the convergence to the H^-noim of (iV, U, 0) at the 
price of losing the decay rate. To be precise, we take the inner product of the system (14.6|) 
with A'^{N, U, 0) to obtain an energy inequality. Among the terms, we consider the typical 
term 

f ft2div((l + p^)Au^) ^ 2 ^ + pf^)Au^) 

J 36(l + p^) 36(l + p^) 

< J|VA0f + ;i4(||VdivAu'‘||i. + IIAp'^lliellAu'^llia 

+ llVp'^lli.lldivAn^lli. + ||Vp'‘||i.||Au'‘||i.) 

< J|VA0f + WdivAu'^lli. + (1 + ||p'^|||,3)||u"||?,3) 

<^||vA0f+ c;i2^ 

thanks to the estimate in (14.41) . The other terms, either depending linearly on {N, [/, 0) or 
depending on the Planck constant h, can be treated similarly by making use of (14.41) . Finally, 
we arrive at the inequality 

[/, emWrn + 0)11?,. < cilKiv, t/, 0)(t)||?,3 + C 2 h^, 

which implies, thanks to the Gronwall inequality, that 

||(iV,t/,0)(t)||?,3 < [c2e^^7ci] 

In particular, ci and C 2 are independent of h. This completes the proof. 
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